/$£# 77ff- /d / ' 3 ''7 


3 1176 01328 6126 ' 


NASA Technical Memorandum 101319 


NAS A-TM- 101319 19880017790 


A Computational Analysis of 
Under-Expanded Jets in the 
Hypersonic Regime 


Andrew T. Hsu 
Sverdrup Technology, Inc. 

NASA Lewis Research Center Group 
Cleveland, Ohio 

and 

Meng-Sing Liou 

National Aeronautics and Space Administration 
Lewis Research Center 
: Cleveland , Ohio 


Presented at the 

Atmospheric Flight Mechanics Conference 

sponsored by the American Institute of Aeronautics and Astronautics 
Minneapolis, Minnesota, August 15-17, 1988 



NOV 3 0 1968 


LIBRARY NASA 
HAMPTON, VIRGINIA 





A COMPUTATIONAL ANALYSIS OF UNDER-EXPANDED 
JETS IN THE HYPERSONIC REGIME 

Andrew T. Hsu* 

Sverdrup Technology, Inc. 

NASA Lewis Research Center Group, Cleveland, Ohio 44135 

and 

Meng-Sing Liou** 

National Aeronautics and Space Administration 
Lewis Research Center, Cleveland, Ohio 44135 


Abstract 

Underexpanded axisymmetric jets are studied nu- 
merically using a full Navier-Stokes solver. Emphasis has 
been given to supersonic and hypersonic jets in supersonic 
and hypersonic ambient flows, a phenomenon previously 
overlooked for the most part. The present work 
demonstrates that the shear layers and shock patterns 
in a jet plume can be captured without complicated vis- 
cous/inviscid and subsonic/supersonic coupling schemes. 
In addition, a supersonic pressure relief effect has been 
identified for underexpanded jets in supersonic ambient 
flows. While it is well known that an underexpanded jet 
in a quiescent ambience (or subsonic ambience) contains 
multiple shock cells, the present study shows that because 
of the supersonic pressure relief effect, an underexpanded 
jet in a supersonic or hypersonic ambience contains only 
one major shock cell. 

I. Introduction 

The phenomenon of under-expanded jets can be 
found in many engineering applications, such as the air- 
craft plume/afterbody interaction, 1 rocket plume signa- 
ture prediction, 2 and jet noise prediction. 3,4 Many exper- 
iments have been carried out to study this problem. 3 ’ 5 ” 7 
Most of these experiments are restricted to underex- 
panded jets in quiescent ambient flows, with the exception 
of Reid and Hastings’s work, 7 which reported experimen- 
tal studies of underexpanded jets in a Mach 2 freestream. 
Chuech et al. 6 studied turbulence structures in underex- 
panded jets. Seiner and Norum 3 reported detailed ex- 
perimental studies on underexpanded jets in a quiescent 
environment; their data is used to calibrate the present 
study. 

Because jet plumes contain shocks, expansion waves, 
and both supersonic and subsonic zones, numerical anal- 
yses of underexpanded jets often involve complicated vis- 
cous/inviscid and subsonic/supersonic coupling schemes. 
Over the years, many numerical schemes have been de- 
vised to capture the intricate shock patterns in underex- 
panded jets. Salas 8 calculated inviscid jet plumes using 
a shock fitting method. Dash et al. 2,0 “ n made extensive 
studies of underexpanded jets using an Euler solver and 
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a parabolized Navier-Stokes solver; their viscous flow re- 
sults compared quite well with Seiner’s experimental data. 
However, these solvers either have to divide the flowfield 
into domains of different flow properties or have to rely 
on viscous/ inviscid iterations, and thus their applications 
are not straightforward. 

In the present work, a full Navier-Stokes equation 
solver is used to study underexpanded jets. The advan- 
tage of the present scheme lies in its simplicity: no vis- 
cous/inviscid iteration is needed, and the subsonic and 
supersonic zones are treated uniformly. Moreover, this 
scheme is capable of capturing the details of the barrel 
shocks and Mach disks in the jet plume. 

To date, most of the numerical studies of underex- 
panded jets are restricted to the low supersonic regime, 
and emphasis has been on jet flows in quiescent environ- 
ments. Salas 8 and Dash 10 made some inviscid flow calcu- 
lations on underexpanded jets in supersonic freestreams, 
but because of the limitations of the models used, the in- 
teraction between the jets and the ambient flows is not 
thoroughly discussed. Since the present scheme treats 
the jet and the ambient flow as one flowfield, the interac- 
tion between the jet and the supersonic ambient flow is 
studied in detail. It is found that this interaction consid- 
erably alters the jet structure. Because it is difficult to 
collect experimental data, little is known about jet flows 
at hypersonic conditions. Thus the present work stud- 
ies the characteristics of hypersonic jet flows numerically, 
and obtains information that is difficult to produce ex- 
perimentally. ’ 

The present work studies only the near jet flowfield 
immediately after the nozzle, where the effects of turbu- 
lence are not predominant. Therefore, only laminar jet 
flow is studied. Finally, cold flow is considered in this 
work; no attempt is made to assess the real gas effect. 

II. Mathematical and Numerical Formulations 

2.1 Equations and solver 

The full Navier-Stokes equations in conservation law 
form are used in the present study: 

d t q + d z E + d y F = Rc~ l (d z R + d v S) 
where q is a vector containing the conservation variables: 
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The vectors E and F are the inviscid flux vectors: 
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generated may not be smooth enough to be used in actual 
calculations. A filter type of smoothing scheme is devised 
to filter out the high frequency oscillations in the grid. 

The present numerical results show that, with the 
same number of grid points, the use of an adaptive grid 
can improve the accuracy of the solution by more than 
ten percent. 


and R and 5 are the viscous flux vectors: 
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An existing computer code, PARC, originally devel- 
oped by Pulliam et al. 12 as ARC, and later modified by 
the AEDC Group of Sverdrup Technology, 13 is used in the 
present study. This code solves the above equations using 
the approximate factorization algorithm originally devel- 
oped by Beam and Warming. Central differencing is used 
to discretize the spacial derivatives, and backward differ- 
encing is used for the time derivatives. To avoid having to 
solve a block pentadiagonal matrix, the Jacobian matrices 
are diagonalized using their eigenvalues and eigenvectors. 
This procedure results in a set of scalar pentadiagonal 
equations. Second-order and fourth-order artificial dis- 
sipation terms are used in the code to ensure stability 
and convergence. The axisymmetric version of the code 
is used in this study. 


III. Numerical Analysis of Jet Flows 

An underexpanded axisymmetric jet contains a fairly 
complicated structure of shocks and expansion waves, and 
Mach disks (normal shocks) can create pockets of sub- 
sonic flows in a supersonic jet. Moreover, the ambient 
flow may be supersonic or subsonic, adding more com- 
plexity to the problem. To accommodate the diversity of 
the flowfield, Dash et al. 0 ” 11 have used a rather sophis- 
ticated numerical scheme. However, as will be shown in 
the present work, the intricacy of an underexpanded jet 
can be treated using a uniform scheme, without having 
to differentiate various flow regions. 

Several different jet Mach numbers, freestream Mach 
numbers, and pressure ratios have been considered in the 
present study. These flow conditions are given in Table 
1. 


case 

Mj 

M a 

Pi/Po 

1 

2 

0 

1.45 

2 



2 

3 

3 

2 

5 

4 



10 

5 



3 

6 

6 

5 

5 

7 



10 


Table 1. Flow conditions 


2.2 Grid 

In order to improve accuracy, an adaptive grid pack- 
age is developed by the first author to accommodate the 
intricate shock patterns and the wave-like shear layer in 
the jet flowfield. 

This grid adaptation package uses an arc equi- 
distribution method suggested by D.A. Anderson. 14 In 
this scheme, the grid size is required to be inversely pro- 
portional to the gradient of flow variables, e.g., in the 
y-direction: 



where Q can be any of the physical quantities such as the 
velocity, pressure, etc. A slight variation from the scheme 
originally proposed by Anderson 14 yields 

Ay,[l + p{dQldy)}\ = const . 

The formula above can be applied to any coordinates to 
generate a flow-adaptive grid. However, the grid thus 


3.1 A Mach 2 jet in a quiescent ambience 

It is known that compressible flow solvers usually 
have difficulties in dealing with very low Mach number 
flows. According to Chuech, 6 a coflow must be added to 
stabilize the solution when using Dash’s SCIPVIS code 
to solve jet flows in a quiescent environment. Our expe- 
rience shows that changes in coflow can change the cell 
length in the jet considerably. In the present study, a 
jet flow into a true quiescent environment is studied. Be- 
cause regions of zero Mach number exist, cpnvergence is 
slow, and is sensitive to the grid and to the boundary 
conditions. Nonetheless, converged solutions can be ob- 
tained for the relatively low jet-to- ambient pressure ratio 
Pj/Pa = 1-45, with a jet Mach number Mj = 2. A typical 
flow-adapted grid used in the present study b shown in 
Figure 1. The use of flow-adaptive grids enables one to 
capture the shocks and the shear layer better, thus im- 
proving the accuracy of the predicted pressure, etc. De- 
tails of the effects of adaptive grids on flow solutions will 
be reported in a separate paper. 

Dash and Wolf 10 gave a thorough description of the 
flow features for underexpanded jets in a quiescent am- 
bience. Figure 2 is the computed Mach number contour 
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from the present study. All the basic features described 
by Dash and Wolf are successfully captured in the present 
numerical results. The calculated jet centerline pressure 
is compared with Seiner’s experimental data 3 in Figure 
3. Except for an under-prediction after the first Mach 
disk, the numerical result agrees fairly well with the ex- 
perimental data. It is worth noticing that the under- 
prediction in the present result is not any worse than 
other numerical results obtained from more sophisticated 
solution schemes. 

3.2 A Mach 3 jet in a Mach 2 freestream 

Since underexpanded jets find most of their applica- 
tions in supersonic flights, the study of jets in supersonic 
ambient flows is probably more important than the study 
ofjets in still air. In the present paper, a Mach 3 under- 
expanded jet in a Mach 2 freestream is studied. Three 
jet-to-ambience pressure ratios, Pj/p a = 2,5, 10, are con- 
sidered. 

The Mach number contours for the case Pj/p a = 10 
are given in Figure 4. One can see that the jet structure 
and the ambient flow pattern are considerably different 
from those of the case with an ambience of still air. When 
the ambient flow is subsonic (including quiescent), the jet 
boundary is that of a constant pressure. All the waves 
are reflected on this boundary and become waves of the 
opposite sense (e.g., compression waves are reflected as 
expansion waves). This is the mechanism that causes os- 
cillations and a series of shock cells in the underexpanded 
jet. However, when the ambient flow is supersonic, the jet 
boundary is no longer a boundary of constant pressure. 
Because the jet expands and the freestream changes direc- 
tion, an oblique shock develops in the ambient flowfield at 
the jet exit, which increases the pressure at the jet bound- 
ary and essentially decreases the effective jet-to-ambience 
pressure ratio. The pressure build-up is relieved when 
the turning jet boundary induces a series of expansion 
waves, which can also be interpreted as waves in the jet 
passing through the jet boundary. When the shock from 
the Mach disk meets the jet boundary, it helps the jet 
flow adjust to the direction of the freestream. There is 
an oblique shock in the freestream originating from the 
point where the jet boundary turns to the freestream di- 
rection. Again, this shock can be explained either as aris- 
ing from the turning of the jet boundary or as a shock 
from the jet passing through the boundary. The above 
described process helps to adjust the jet pressure to that 
of the ambience very quickly, normally within one shock 
cell, although the cell lengths are much longer here than 
in the still air cases. This mechanism of adjusting the 
jet pressure to the ambient pressure is what we call the 
supersonic pressure relief effect. 

The jet boundary pressure for Pj/pa = 5 is given 
in Figure 5, and the jet centerline pressures for all three 
supersonic flow cases are given in Figure 6. The centerline 
pressure distributions show that, in all the three cases, 
the pressure oscillations are negligible after the first shock 
cell. Dash and Wolf 10 calculated a Mach 2 jet in a Mach 
2 freestream; similar phenomenon was observed, but no 
explanations were given in their study. 


3.3 A Mach 6 jet in a Mach 5 freestream 

As the title of this paper indicates, the primary goal 
of the present work is to understand the structures of 
hypersonic jets. Towards this end, a Mach 6 jet in a Mach 
5 freestream is studied. Three jet-to-ambience pressure 
ratios are considered; they are pj/p a = 3, 5, and 10. 

The present results show that the basic structure of 
a hypersonic jet is similar to that of a supersonic jet. 
(Here ‘supersonic* and ‘subsonic* also refer to the ambient 
flows.) All the previous analysis about supersonic jets also 
applies here. The major differences for the two different 
flight regimes are that, in the hypersonic cases, the cell 
lengths are much longer, and the shocks and expansion 
waves induced in the freestream by the jet stay fairly close 
to the boundary of the jet. A Mach number contour for 
the hypersonic jet with an underexpansion pressure ratio 
of pj/pa = 5 is shown in Figure 7. 

The centerline pressure distributions for the three hy- 
personic jet cases mentioned earlier are given in Figure 
8. A comparison between Figures 6 and 8 shows that as 
the flight Mach number increases, the oscillation in the 
centerline pressure decreases. This result occurs because 
a higher Mach number causes a stronger shock at the 
jet exit, and thus the supersonic pressure relief effect is 
stronger there. 

The cell lengths of the underexpanded jets as a func- 
tion of the jet-to-freestream pressure ratio are given in 
Figure 9, where the present numerical data are being com- 
pared with the experimental data given by Love, et al., 15 
and Chiang. 5 Here the relative shock cell length is defined 
as the distance between the nozzle exit and the first Mach 
disk, /*<*, divided be the nozzle exit radius 72/. The ex- 
perimental data are for underexpanded jets in quiescent 
ambiences. Figure 9 shows that with the same relative ve- 
locity between the jet and the freestream, the cell lengths 
can be significantly different for various flight Mach num- 
bers. 

IV. Concluding Remarks 

The present numerical study demonstrates that the 
phenomenon of underexpanded jets can be analyzed us- 
ing full Navier-Stokes solvers without resorting to special 
treatments of the jet flowfield. It also shows that the jet 
structures are considerably different for subsonic and su- 
personic ambient flows. A supersonic pressure relief effect 
has been identified. This pressure relief effect reduces the 
pressure oscillation in an underexpanded jet. The study 
also shows that the length scales in hypersonic jets are 
considerably different from those in supersonic jets. 
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Figure 1. Flow adapt grid for an underexpanded jet. 



Figure 2. Mach number contour for an unerexpanded jet. 
My = 2, M q = 0, Py / p a = 1.45. 


4 







Figure 4. Mach number contour for a Mach 3 jet in a 
Mach 2 ambient flow. pj/p a = 10. 



Figure 5. Jet boundary pressure distribution for a Mach 
3 jet in a Mach 2 freestream. py/Po = 5- 
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Figure 6. Jet centerline pressure distributions for Mach 3 
jets in mach 2 freestreams. py/p 0 = 2, 5, and 10. 



Figure 7. The Mach number contour for a hypersonic 
underexpanded jet in a hypersonic freestream. My = 6, 
M fl = 5, p 3 /pa = 5. 


6 








MSA Report Documentation Page 

Space Administration 


1. Report No. 2. Government Accession No. 

NASA TM— 1 01319 


4. Title and Subtitle 

A Computational Analysis of Under-Expanded Jets 
In the Hypersonic Regime 


3. Recipient's Catalog No. 



6. Performing Organization Code 


7. Authors) 


Andrew T. Hsu and Meng-Sing Llou 


8. Performing Organization Report No. 

E— 43 1 7 

10. Work Unit No. 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


505-62-81 

11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Presented at the Atmospheric Flight Mechanics Conference sponsored by the 
American Institute of Aeronautics and Astronautics, Minneapolis, Minnesota, 
August 15-17, 1988. Andrew T. Hsu, Sverdrup Technology, Inc., NASA Lewis 
Research Center Group, Cleveland, Ohio 44135; Meng-Sing Llou, NASA Lewis 
Reseats 


16. Abstract 

Underexpanded axisymmetric jets are studies numerically using a full Navier- 
Stokes solver. Emphasis has been given to supersonic and hypersonic jets In 
supersonic and hypersonic ambient flows, a phenomenon previously overlooked for 
the most part. The present work demonstrates that the shear layers and shock 
patterns In a jet plume can be captured without complicated viscous/lnviscid and 
subsonic/supersonic coupling schemes. In addition, a supersonic pressure relief 
effect has been identified for underexpanded jets In supersonic ambient flows. 
While It Is well known that an underexpanded jet In a quiescent ambience (or 
subsonic ambience) contains multiple shock cells, the present study shows that 
because of the supersonic pressure relief effect, an underexpanded jet in a 
supersonic or hypersonic ambience contains only one major shock cell. 


17. Key Words (Suggested by Authors)) 

Jet plume - 
Numerical analysis 
Hypersonlcs 


18. Distribution Statement 


Unclassified - Unlimited 
Subject Category 02 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No of pages 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 











FOURTH CLASS MAIL 


National Aeronautics and 
Space Administration 

Lewis Research Center 

Cleveland, Ohio 44135 


ADDRESS CORRECTION REQUESTED 




Official Business 
Penalty for Private Use $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA 451 




